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I.  INTRODUCTION 


Contract  N00014-K-0574,  which  began  1  August  1984,  provides  for 
basic  research  on  nonlinear  acoustics.  The  contract  is  the  successor  to 
Contract  N00014-75-C-0867  (ended  31  August  1984)  and  also  continues 
research  begun  under  Contract  N00014-82-K-0805  (ended  30  November 
1984)  on  the  specific  topic  of  nonlinear  effects  in  long  range  underwater 
propagation.  The  period  of  this  report  is  1  August  1984  -  31  October  1985. 
However,  because  much  of  the  current  work  began  under  the  previous  two 
contracts,  for  clarity  this  report  covers  some  of  the  work  done  prior  to  the 
present  contract.  The  basic  objective  of  the  research  is  to  investigate  the 
behavior  of  finite-amplitude  sound,  in  particular  to  investigate  effects  that 
do  not  occur  in  small -signal  sound  fields. 

The  personnel  involved  in  the  research  are  as  follows. 

Senior  personnel 

D.  T.  Blackstock,  principal  investigator 

R.  Buckley,  consultant.  Institute  of  Sound  and  Vibration  Research 
(ISVR),  University  of  Southampton,  UK. 

C.  L.  Morfey,  consultant,  ISVR,  University  of  Southampton,  UK. 

W.  M.  Wright,  consultant,  Physics  Dept.,  Kalamazoo  College, 
Kalamazoo,  MI. 

Postdoctoral  fellow 

M.  F.  Hamilton,  Ph.D.,  Pennsylvania  State  University  (1983). 

Graduate  students 

F.  D.  Cotaras,  M.S.  student  in  Electrical  Engineering.  Completed 
M.S.  degree  August  1985;  Ph.D.  student  beginning  September 
1985. 

D.  A.  Nelson,  M.S.  student  in  Mechanical  Engineering.  Completed 

M.S.  degree  December  1985. 

J.  A.  TenCate,  Ph.D.  student  in  Mechanical  Engineering. 


II.  PROJECTS 


1.  Nonlinear  effects  in  long  range  underwater  propagation  (Cotar as, 
Morfey,  and  Buckley).  This  work  began  under  Contract  N00014-82-K-0805. 
The  research  was  started  by  C.  L.  Morfey,  who  was  at  ARLsUT  on  leave 
from  ISVR  during  the  first  six  months  of  1983.  Assisting  as  a  consultant  was 
Richard  Buckley,  also  of  ISVR,  who  functioned  primarily  under  the  direction 
of  Morfey.  After  Morfey  returned  to  ISVR  in  July  1983,  he  continued  the 
research,  at  a  much  reduced  level  of  effort,  as  a  consultant  to  ARLtUT.  In 
September  1983,  F.  D.  Cotaras,  at  ARLsUT  on  educational  leave  from  DREA 
in  Canada,  began  work  on  the  project  for  a  master's  degree.  A  brief 
description  of  the  work  done  prior  to  1  August  1984  will  be  followed  by  a 
report  of  progress  since  that  date. 

Scientists  and  engineers  working  on  problems  in  long  range  under¬ 
water  propagation  depend  almost  entirely  on  linear  acoustical  theory.  They 
recognize  that  in  order  to  produce  signals  that  travel  great  distances  the 
source  may  have  to  generate  a  sound  of  finite  amplitude.  They  assume, 
however,  that  after  the  signal  has  propagated  a  certain  range,  say  one 
hundred  or  a  few  hundred  meters,  the  amplitude  will  have  decreased  to  the 
point  that  nonlinear  effects  are  no  longer  important.  Our  research  was 
guided  by  two  questions.  (1)  Is  there  really  a  range  beyond  which 

nonlinearity  may  be  ignored?  (Certain  evidence  pointed  to  a  negative 
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answer  to  this  question.  ’  )  (2)  If  such  a  range  does  exist,  what  is  its 

magnitude  and  on  what  does  it  depend? 


H.  W.  Marsh,  R.  H.  Mellen,  and  W.  L.  Konrad,  "Anomalous  absorption  of 
pressure  waves  from  explosions  in  sea  water,"  J.  Acoust.  Soc.  Am.  38, 
326-328  (1965). 

D.  A.  Webster  and  D.  T.  Blackstock,  "Experimental  investigation  of  outdoor 
propagation  of  finite-amplitude  noise,"  NASA  Contractor  Report  2992, 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin  (August 
1978)  (N78-31876). 


The  three  linear  theories  most  commonly  used  in  long  range  propaga¬ 
tion  studies  are  ray  theory,  the  parabolic  approximation,  and  normal  mode 
theory.  For  three  reasons  we  selected  ray  theory  for  our  work.  First,  ray 
theory  is  simple  and  versatile.  Second,  much  was  known  about  how  to 
extend  ray  theory  to  cover  finite-amplitude  waves.  Third,  the  prospects  for 
analytical  solutions  in  certain  cases  were  good. 

When  Morfey  began  the  research,  he  divided  it  into  three  tasks. 

I.  Shock  pulse  propagation  in  a  homogeneous  ocean.  The  intent 
here  was  to  develop  a  computer  program  to  calculate  the 
propagation  of  sound  from  realistic  explosions  in  a  homo¬ 
geneous  but  otherwise  realistic  ocean.  In  other  words,  we 
wished  to  implement  weak  shock  theory  for  spherical  waves  in 
a  relaxing,  viscous  fluid. 

II.  Nonlinear  propagation  in  a  depth  dependent  ocean.  Shocks  and 
dissipation  were  ignored  in  this  task.  Of  primary  interest  was 
the  effect  of  inhomogeneity  on  nonlinear  distortion.  An 
analytical  model  based  on  nonlinear  geometrical  acoustics  was 
to  be  developed. 

III.  Nonlinear  propagation  in  a  caustic  region.  This  topic,  which 
required  special  attention  because  of  the  confluence  of  three 
effects  —  nonlinear  distortion,  focusing,  and  diffraction  —  was 
assigned  to  consultant  Richard  Buckley. 

The  overall  plan  was  to  work  on  each  of  the  three  tasks  separately  and, 
when  they  were  completed,  to  put  the  results  all  together,  probably  in  the 
form  of  a  single  computer  program.  One  would  then  be  able  to  calculate  the 
propagation  of  an  intense  underwater  signal  over  great  distances,  under  a 
variety  of  conditions. 


Starting  from  where  Morfey  left  off  on  Tasks  I  and  II,  Cotaras  carried 
out  further  work  on  the  two  tasks  as  needed  and  combined  the  results.  In 
retrospect  this  was  not  a  small  order,  despite  the  fact  that  Morfey  had  in 
one  sense  nearly  completed  Task  n.^ 

During  the  period  covered  by  this  report  Cotaras  accomplished  four 
things.  First,  he  completed  all  the  analytical  work  to  show  how  the 
equations  of  nonlinear  geometrical  acoustics  are  derived  from  the  general 
equations  of  fluid  mechanics  for  an  inhomogeneous  medium.  Second,  he 
completed  the  computer  program  (which  had  been  started  by  Morfey)  by 
which  propagation  calculations  based  on  nonlinear  geometrical  acoustics  are 
carried  out.  The  program  represents  a  merger  of  Tasks  I  and  II.  He  tested 
the  program  by  running  it  for  cases  for  which  analytical  results  are  known. 
Third,  he  applied  the  program  to  two  different  source  signals  (explosions)  of 
interest.  A  parametric  study  of  propagation  up  to  57  km  was  performed. 
An  interim  report  was  given  at  the  Fall  1984  Acoustical  Society  of  America 
(ASA)  Meeting  in  Minneapolis  (84-4).^  One  of  the  unique  features  of  the 
program  is  the  option  available  to  the  operator  of  turning  off  nonlinear 
effects  beyond  a  given  range.  This  option  was  very  useful  in  helping  to 
answer  the  question  with  which  we  began  this  project:  Beyond  what  range  is 
it  safe  to  ignore  nonlinear  effects?  The  answer  is  not  simple,  but  we  did 
learn  that  the  "cutoff  range"  increases  with  the  frequency  band  in  which  one 
is  interested.  Finally,  Cotaras  wrote  his  thesis  (85-7),  the  abstract  of  which 
is  quoted  here. 

"In  this  thesis  the  propagation  of  finite  amplitude  acoustic  signals 
through  an  inhomogeneous  ocean  is  investigated  both  analytically  and 
numerically.  The  effects  of  reflections  and  focusing  are  not  considered. 
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C.  L.  Morfey  "Nonlinear  Propagation  in  a  Depth -Dependent  Ocean,"  Tech¬ 
nical  Report  ARL-TR-84-11,  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  1  May  1984  (ADA  145  079). 


Numbers  given  in  this  style  refer  to  items  in  the  Chronological  Bibli¬ 
ography,  e.g.,  84-4  means  the  fourth  entry  in  the  list  for  1984. 


From  simplified  versions  of  the  lossless  hydrodynamics  equations  the 
theories  of  linear  and  nonlinear  geometrical  acoustics  are  developed.  Losses 
are  accounted  for  directly  in  the  numerical  routine.  The  eikonal  equation, 
from  which  an  equation  for  the  ray  paths  is  derived,  is  assumed  to  be  the 
same  for  both  small-signal  and  finite  amplitude  waves.  The  transport 
equation  is  found  to  be  different,  however.  The  transport  equation  leads  to 
a  standard  first-order  progressive  wave  equation,  linear  for  small-signal 
waves,  but  nonlinear  for  finite  amplitude  waves.  All  the  analysis  is  carried 
out  in  the  time  domain  and  is  for  a  fully  inhomogeneous  ocean. 

"In  the  numerical  study  the  ocean  is  assumed  to  be  stratified.  The 
effects  of  inhomogeneity,  ordinary  attenuation  and  dispersion,  and  nonlinear 
propagation  are  investigated  using  a  numerical  implementation  of  nonlinear 
geometrical  acoustics.  Two  explosion  waveforms  are  considered:  a  weak 
shock  with  an  exponentially  decaying  tail  and  a  more  realistic  waveform  that 
includes  the  first  bubble  pulse.  Numerical  propagation  of  the  simpler  wave 
along  a  58.1  km  path  starting  at  a  depth  of  300  m  leads  to  the  following 
conclusions:  (1)  The  effect  of  inhomogeneity  on  nonlinear  distortion  is 
small.  (2)  Dispersion  plays  an  important  role  in  determining  the  arrival  time 
of  the  pulse.  (3)  Neither  nonlinearity  nor  ordinary  attenuation  (and  disper¬ 
sion)  are  paramount;  both  need  to  be  included.  For  the  more  realistic  wave 
the  propagation  is  along  a  23  km  ray  path  starting  from  a  depth  of  4300  m. 
Two  charge  weights,  0.818  kg  and  22.7  kg  TNT,  are  used.  In  each  case  the 
energy  spectrum  of  the  signal  obtained  by  considering  finite  amplitude 
effects  for  the  entire  23  km  path  is  compared  with  spectra  obtained  by 
neglecting  finite  amplitude  effects  (1)  entirely,  (2)  after  the  first  150  m,  and 
(3)  after  the  first  1100  m.  Finite  amplitude  effects  are  found  to  be  of  small 
consequence  in  the  case  of  the  0.818  kg  TNT  explosion  for  frequencies 
below  6  kHz  at  distances  beyond  1100  m.  For  the  22.7  kg  explosion  the 
corresponding  quantities  are  4  kHz  and  1100  m." 

Morfey  and,  to  a  lesser  extent,  Buckley  have  continued  their  work  on 
the  project  as  consultants.  Most  of  the  effort  has  been  toward  writing 
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journal  articles  on  work  that  has  been  accomplished.  In  November  1984 
Buckley  submitted  a  manuscript  on  his  work  on  caustics  to  the  Journal  of 
the  Acoustical  Society  of  America.  As  a  result  of  comments  by  reviewers, 
he  now  plans  to  divide  the  manuscript  into  two  articles,  one  on  the 
time-domain,  linear  theory  of  caustics  and  the  other  on  nonlinear  theory.  A 
paper  on  the  latter  was  given  at  the  Spring  1985  ASA  Meeting  in  Austin 
(85-3).  Morfey  is  also  working  on  or  has  plans  for  two  journal  articles.  One 
article  is  based  on  his  Task  I  results,  which  were  reported  partly  at  the  Kobe 
International  Symposium  on  Nonlinear  Acoustics  (ISNA)^  and  more  recently 
at  the  Austin  ASA  Meeting  (85-2).  The  second  article  will  be  based  on  his 
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technical  report. 

2.  Nonlinear  effects  in  air-filled,  bulk  porous  materials  (Nelson).  This 


work  began  in  1978  under  NASA  Grant  NSG  3198.  The  first  phase  was 
completed  in  1982,^  at  which  time  Nelson  began  working  on  the  problem. 
The  project  continued  to  receive  its  main  support  from  the  NASA  grant  until 
its  expiration  28  February  1985.  Nelson  completed  his  thesis  in  December 
1984  (84-9),  which,  with  an  added  chapter  on  lined  ducts,  became  a 
technical  report  (85-1).  The  following  is  the  abstract  from  the  report. 

"The  subject  of  this  investigation  is  the  propagation  of  high  intensity 
sound  waves  through  an  air-filled  porous  material.  The  material  is 
assumed  (1)  to  be  rigid,  incompressible,  and  homogeneous,  and  (2)  to  be 
adequately  described  by  two  properties:  resistivity  r  and  porosity  12  .  The 


C.  L.  Morfey,  "Aperiodic  signal  propagation  at  finite  amplitude:  some 
practical  applications,"  10th  International  Symposium  on  Nonlinear 
Acoustics,  Kobe,  Japan,  24-28  July  1984. 


H.  L.  Kuntz,  "High-intensity  sound  in  air  saturated  fibrous  bulk  porous 
materials,"  Ph.D.  dissertation,  The  University  of  Texas  at  Austin.  Also 
issued  under  the  same  name  as  Technical  Report  ARL-TR-82-54,  Applied 
Research  Laboratories,  The  University  of  Texas  at  Austin  (3  September 
1982)  (ADA  121  450). 


resistivity  was  measured  as  a  function  of  velocity  for  static  flows  and  found 
to  follow  the  empirical  relation  r^+^usgnfu).  This  relation  is  assumed  to 
apply  for  acoustic  signals  as  well.  Ordinary  hydrodynamic  nonlinearity 
(which  leads  to  shock  formation)  is  neglected  because  of  the  very  high 
attenuation  in  the  porous  material.  The  resulting  wave  equation  is  still 
nonlinear,  however,  because  of  the  u  sgn(u)  term  in  the  resistivity.  The 
equation  is  solved  in  the  frequency  domain  as  an  infinite  set  of  coupled 
inhomogeneous  Helmholtz  equations,  one  for  each  harmonic.  An  approxi¬ 
mate  first  integral  formulation  of  these  equations  gives  relations  for 
progressive  waves.  The  source  wave  considered  is  a  slightly  distorted 
intense  tone,  that  is,  a  finite-amplitude  fundamental  accompanied  by  weak 
higher  harmonics.  An  approximate  but  analytical  solution  leads  to  predic¬ 
tions  of  excess  attenuation,  saturation,  and  phase  speed  reduction  for  the 
fundamental  component.  A  more  general  numerical  solution  is  used  to 
calculate  the  propagation  curves  for  the  higher  harmonics.  The  u  sgn(u) 
nonlinearity  produces  a  cubic  distortion  pattern;  when  the  input  signal  is  a 
pure  tone,  only  odd  harmonic  distortion  products  are  generated.  Quantita- 
tive  experiments  were  performed  on  batted  Kevlar  w  29  having  porosities  in 
the  range  fi=0.94  to  0.98.  Measurements  were  taken  in  the  400  to  6200  Hz 
frequency  range  for  small-signal  waves,  500-1500  Hz  for  finite-amplitude 
waves.  Source  levels  for  the  finite-amplitude  waves  extended  to  160  dB. 
Qualitative  confirmation  of  predictions  has  been  obtained  in  all  cases,  and 
quantitative  confirmation  in  most  cases.  The  theory  is  then  applied  to  the 
problem  of  propagation  and  attenuation  in  lined  ducts.  A  vector  version  of 
the  wave  equation  is  derived  along  with  the  corresponding  inhomogeneous 
Helmholtz  equations.  These  are  used  as  the  basis  of  a  perturbation  solution 
for  reflection  from  a  porous  half -space  and  propagation  in  lined  ducts.  The 
solution  becomes  extremely  complicated,  however,  when  used  for  high 
intensity  waves  in  lined  ducts.  An  ad  hoc  model  is  therefore  presented  for 
the  purpose  of  illustrating  the  gross  nonlinear  effects  on  absorption." 

A  report  on  the  later  stages  of  the  work  was  given  at  the  Spring  1985 
ASA  Meeting  in  Austin  (85-4). 


3.  Nonlinear  interaction  in  a  rectangular  waveguide.  (Hamilton  and 


TenCate).  This  work  is  an  outgrowth  of  Hamilton's  Ph.D.  research  (84-6,7) 
and  TenCate's  M.S.  research  (84-2).  Experiments  and  theoretical  calcula¬ 
tions  were  carried  out  for  two  different  excitations  of  the  waveguide, 
referred  to  here  for  convenience  as  Cases  A  and  B.  In  Case  A  a  source 
launches  a  monochromatic  wave  called  the  fundamental  (of  frequency  f) 
that  travels  down  the  waveguide  in  the  1,0  mode.  The  discussion  that 
follows  may  most  easily  be  described  in  terms  of  small-signal  concepts.  In 
the  1,0  mode  the  wave  bounces  back  and  forth  between  the  top  and  bottom 
walls  of  the  waveguide  as  shown  in  Fig.  1  (the  0,0  mode  is  also  shown;  it 
plays  no  role  in  Case  A).  The  angle  the  wave  makes  with  the  waveguide 
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Fig.  1.  Propagation  in  the  0,0  and  1,0  modes  of  a  rectangular  waveguide 
(f  =  frequency,  cq  =  small-signal  sound  speed,  a  =  height  of  the 
guide). 


axis  is  given  by  the  formula  in  the  figure.  As  the  frequency  is  lowered 

O 

toward  cutoff,  6  approaches  90  .  Far  above  cutoff  on  the  other  hand,  6 

O 

approaches  0  .  The  bouncing  wave  may  be  decomposed  into  upward  and 
downward  traveling  waves  that  make  angles  +  6  and  -  e,  respectively,  with 
the  waveguide  axis.  Now,  when  the  wave  is  of  finite  amplitude,  it  generates 
a  family  of  harmonic  distortion  components.  Resonant  or  synchronous 


interaction  (like  that  which  takes  place  when  a  plane  wave  propagates  in  an 
open  medium)  occurs  when  all  the  harmonics  in  the  family  travel  at  the 
same  speed  and  in  the  same  direction.  This  happens  when  the  second 
harmonic  travels  in  the  2,0  mode,  the  third  harmonic  in  the  3,0  mode,  and  so 
on.  The  angle  d  is  then  the  same  for  all  members  of  the  family  (all 
members  have  the  same  group  velocity).  It  turns  out  that  the  analysis  is 
simplest  when  the  frequency  is  such  that  6  is  in  the  middle  range,  that  is, 
neither  close  to  cutoff  nor  close  to  high  frequency  operation.  In  this  range 
the  upward  and  downward  traveling  waves  do  not  have  much  effect  on  each 
other.  Near  cutoff  (0  —  90°)  the  standing  wave  nature  of  the  field  becomes 
important  and  the  interaction  between  the  upward  and  downward  waves 

O 

cannot  be  ignored.  At  high  frequency  (0  —  0)  the  two  waves  begin  to 
reinforce  each  other  strongly.  Theory  and  data  for  a  middle  range 
experiment  is  shown  in  Fig.  2.  Boundary  layer  losses  at  the  tube  walls  have 
been  accounted  for.  It  is  seen  that  the  theoretical  predictions  are  borne  out 
by  the  measurements.  This  work  was  reported  at  the  Austin  ASA  Meeting 
(85-5). 


Fig.  2.  Computed  and  measured  propagation  curves  for  the  fundamen¬ 
tal  and  second- fifth  harmonic  distortion  components  in  a 


rectangular  waveguide.  The  mode  for  the  fundamental  is  1,0, 
the  frequency  is  3000  Hz,  and  the  angle  0is  55.2°.  The  cutoff 
frequency  for  this  waveguide  is  2464  Hz. 


In  Case  B  two  tones  interact,  one  a  low  frequency  wave  in  the  0,0 
mode,  the  other  a  high  frequency  wave  in  the  1,0  mode  (see  Fig.  1).  The 
distortion  components  of  chief  interest  in  this  case  are  the  sum  and 
difference  frequency  waves.  Dispersion  causes  a  scalloped  pattern  in  the 
amplitude  versus  distance  curves.  The  results,  both  theoretical  and  experi¬ 
mental,  shed  light  on  some  of  the  questions  raised  by  TenCate's  thesis  (84-2), 
for  example,  the  dependence  of  the  coefficient  of  nonlinearity  /?  on  the 
angle  6  at  which  the  two  tones  interact.  Some  of  this  work  was  reported  at 
the  Austin  ASA  Meeting  (85-5).  More  will  be  reported  at  the  Fall  1985  ASA 
Meeting  in  Nashville  (85-8),  and  a  journal  article  is  in  preparation  (85-9). 

4.  Miscellaneous.  A  paper  by  Blackstock  on  a  generalized  Burgers 
equation  was  published  in  the  Journal  of  the  Acoustical  Society  of  America 
(85-6).  This  work  had  been  reported  orally  earlier.7  An  oral  paper  on  a 
topic  in  the  history  of  nonlinear  acoustics  was  given  at  the  Fall  1984 
Meeting  of  the  Texas  Section  of  the  American  Association  of  Physics 
Teachers  in  Houston  (84-8).  Two  earlier  articles  on  nonlinear  acoustics  were 
reprinted  in  Nonlinear  Acoustics  in  Fluids,  the  final  volume  in  the  series 


Benchmark  Papers  in  Acoustics  (84-10  and  84-11).  A  journal  article  by 
Wright  on  excitation  of  a  cylindrical  cavity  by  a  line  thermoacoustic  source 
is  in  preparation  (85-10). 

Two  new  projects  were  just  underway  as  the  report  period  ended, 
both  Ph.D.  dissertation  research  topics.  TenCate  has  begun  work  on  an 
acoustical  chaos  experiment,  intense  standing  waves  in  a  closed  tube.  His 
initial  task  is  to  search  for  subharmonics  in  the  response  of  the  tube  when  it 
is  driven  in  a  higher  order  longitudinal  mode.  For  his  dissertation  topic, 
Cotaras  has  chosen  nonlinear  effects  in  refraction  and  transmission  of  sound 
incident  obliquely  at  an  interface  between  two  media,  for  example,  the 
ocean-sediment  interface.  At  this  point  only  a  literature  survey  has  been 
started. 


D.T.  Blackstock,  "Generalized  Burgers  equation  for  plane  waves,"  9th 
International  Symposium  on  Nonlinear  Acoustics,  20-24  July  1981,  Leeds, 
England. 
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